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Abstract. We derive the amplification of the cosmological magnetic field associated with forming gravitational 
structure. The self-similar solutions of magnetohydrodynamic equations are computed both in linear and nonlinear 
regimes. We find that the relatively fast magnetic field enhancement becomes substantial in the nonlinear phase. 
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1. Introduction 



The hypothesis describing the dynamical role of the pri- 
mordial magnetic field in the formation and evolution of 
gravitational structure frequently occurs in the literature 
(e.g. Wasserman 1978; Kim et al. 1996; Peebles 1995). On 
the contrary, the inverse trend i.e. the amplification of the 
magnetic field during density perturbation collapse is not 
often represented in the context of large-scale structure 
formation. Common practice is to refer to the constraints 
for magnetic field amplification set by the density of col- 
lapsed matter (e.g. Zeldovich et al. 1980). The hints of 
magnetic field existence on cosmological scales excites in- 
terest not only in the absolute value of a frozen field mag- 
nification but also in its growth rate. The early nonlinear 
and previrial phase of the gravitational formation is of par- 
ticular importance, since it results in several megaparsec 
structures observed as superclusters or filaments. An un- 
derstanding of the amplification rate of the structure is 
obviously related to the explanation of the appearance of 
sufficiently strong frozen-in magnetic fields expected at 
this stage of collapse. 

We investigate the mildly nonlinear collapse of cylin- 
drical gravitational structure and give the growth rate of 
the primordial magnetic field as a function of the accretion 
velocity field. The magnetic field growth proceeds inten- 
sively during the phase of fluid compression. The mag- 
netic flux for collapsing plasma is conserved and the mag- 
netic strength changes according to the induction equa- 
tion. It clearly shows that substantial amplification oc- 
curs for strongly compressing flows i.e. for growing divv 
— analogously to the shock processes. The general expres- 
sion for div V is obtained thanks to the self-similar form of 
the hydrodynamic equation. The self-similar presentation 



of magnetohydrodynamics becomes possible in the case of 
rapid density contrast and velocity evolution, when v and 
6 (X a"' [n > 1) i.e. when the Lorentz force neglect in the 
Euler equation is justified. 

The plan of the paper is as follows. In Section 2 we give 
the magnetohydrodynamic (MHD) equations for cylindri- 
cal structures in comoving coordinates of a flat universe. 
We also discuss the assumptions and symmetries allow- 
ing us to separate the induction equation. In Section 3 
we derive the self-similar set of hydrodynamic equations. 
Its linearization and the subsequent comparison with the 
known velocity solutions are also given. We obtain the 
general nonlinear relation divv versus 5. The rate of mag- 
netic field amplification is discussed in Section 4 both in 
the linear and nonlinear regime. In the former we present 
the analytical expression for amplification and the results 
of numerical integration in the latter. 

2. Basic equations and perturbation scheme for 
magnetized fluid 

The magnetic field modifies the standard equations for 
density perturbations and peculiar velocities. These equa- 
tions completed by the Faraday equation for the mag- 
netic field have been extensively studied in the context of 
magnetic field influence on galaxy formation (e.g. Lesch 
& Chiba 1995). The linear approximation of the MHD 
equations has been analyzed in many papers. In the ex- 
treme approach, the primordial magnetic field is respon- 
sible for the origin of density fluctuations (Wasserman 
1978). In his paper the compression of matter is induced 
by the Lorentz force but the linearized hydrodynamic 
equations keep both p and v as small linear values. There 
is no fluid back reaction on the magnetic field. 

In this paper we apply a different approach. The weak 
initial and uniform magnetic field B is embedded in the 
forming large-scale structure which evolves in the dust 
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era. The gravitational collapse changes the strength of the 
magnetic field in the linear and then nonlinear regime. We 
are interested in the nonlinear amplification rate of the or- 
dered field. 

The MEID equations are in comoving coordinates given 
by (e.g. Wassermann 1978) 



|? + 3-V.v + iv-(pv)=0, 
at a a 

dv d 1 , 1„ , 1 ,„ 

^ + -V + -(v-V)v = — - Bx(VxB), 

at a a a Airap 



J-AV' = 47rG(p - po), 

aB a 1 . 

+ 2-B = -Vx vxB , 
at a a 

V • B = 0, 



(1) 

), 
(2) 
(3) 

(4) 
(5) 



where p{t,v) is the matter density {po{t) means the back- 
ground density), v(t, r) is the peculiar velocity of the mat- 
ter forming the structure and tp{t,r) is the density per- 
turbation potential; a{t) is the scale factor given in units 
a{to)- For an assumed fiat FRW cosmological model we 
take the initial value of the uniform magnetic field Bq at 
the initial time to. 

The temporal evolution analyzed here is restricted to 
the after recombination phase and we do not discuss the 
magnetogenesis processes assumed to occur prior to this 
phase. The small value of the initial magnetic energy den- 
sity (relative to the matter density) is assumed. Since the 
field uniformity is maintained within the scale of the form- 
ing structure, it does not perturb the background isotropy. 
The evolving structure and thus the considered field have 
axial symmetry. The velocity of the accreting fiuid is taken 
for simplicity as v = [vr{t,r),vg{t,r),0]. The assumed 
cylindrical symmetry is consistent with the isotropy and 
homogeneity of the FRW model for all length scales large 
with respect to the coherence scale. Finally, the nonlinear 
perturbation equations for density fluctuations of presure- 
less matter, velocity field and the magnetic field B are 
given by 



|J + -V-v+-V-(5v)=0, 
at a a 



(6) 



dv d 1 , 1 „ , 1 /■„ ^ . 
+ -V + -{v ■ V)v = — VVj - (Bo X V X b 

at a a a Anap 



-Kb X V X b), 

J-A?/) = 4ttG{p - po), 

dB „d„ 1„ , 

+ 2-B = -Vx vxB , 
at a a 



(7) 

(8) 

(9) 



where S{t.r) = [p — Po)/Po(t) is the matter density con- 
trast and b(t, r) is the magnetic field perturbation i.e. 
B(t, r) = 'Bo{t)+h{t, ry. In the epoch when the condition 



^ For simplicity, it is assumed here that the background field 
is represented by its uniform component. Following Barrow et 
al. (1997), its sufficiently small initial value (i.e. Bo < 10~^ G) 



of low energy density associated with the magnetic field is 
vahd i.e. B^/p <gC 1, the evolutionary perturbation scheme 
may be simplified. The last expression in the Euler equa- 
tion consists of two terms; the first may be ignored due to 
the small value of Bq and the second is proportional to 
Vp oc V|b|, as can be seen from eqs. (11) and (20) below. 
We shall assume in the first approximation that it is also 
small compared to the density potential gradient. Thus 
in the first iteration step we postulate a force-free mag- 
netic field. The Faraday equation describes the magnetic 
field scaling as a function of accretion velocity. Both the 
amplified magnetic field and the density compression are 
coupled with the velocity field but the induction equation 
decouples from the hydrodynamic equations. Therefore its 
formal solution may be given independently. Applying the 
assumed symmetries and transforming the r.h.s. of eq. (9) 
one obtains the induction equation in the integrable form. 



dBo^ 
dt 



dt 



+ (I-V) (Ba^) = - (Ba^) V- (J 



Its formal solution is 
Ba^ = Bq exp 



dt 



(10) 



(11) 



where v as well as p and V are given by the solutions 
of eqs. (6 8). The successive iterative procedure will cou- 
ple in the next steps the Faraday and gravitational per- 
turbations equations, including the magnetic force to the 
r.h.s. of Euler equation. However this is not essential be- 
fore |B| « ip. 

3. Nonlinear cylindrical perturbation 

The dynamics of collapsing structure modeled from a pres- 

sureless fluid is governed by eqs. (6-8). Applying the as- 
sumed form of accretion velocity, the dynamic equations 
for the fluid transform the continuity equation into 



dS 1 drVr 1 / 6 d(rvr) d6\ 

a7 + ^ + - -^■^+^^^ =0' 

at ar or a \r or or J 

two components of the Euler equation become 



(12) 



dvr d 1 f Vr d rVr 

dt a ^ a \ r dr 



1 



3 (r\2 + r^j) 



1 d^p 
a dr ' 
(13) 

dvg d Ivr drug 

-^ + -ve + ^ = 0. 14 

dt a a r dr 

The peculiar gravitational potential satisfies the Poisson 
equation 

J_d_ / ^\ 
ra? dr \ dr J 

where a{t) = iirGpo. 

does not affect either the background dynamics or the CMB 
spectrum. The obtained final results on B distribution are not 
expected to depend on this (unrealistic on a few Mpc scale) 
assumption. 



(15) 
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Many authors (see for example, Lifshitz 1946; Groth 
& Peebles 1975; Peebles 1980) have studied the analyt- 
ical solutions for the linearized problem of perturbation 
evolution in an expanding Universe. In several cases the 
nonlinear solutions have also been developed (Zeldovich 
1970). Yet, in general, the high degree of complexity of 
the nonlinear set of partial differential equations does not 
allow one to obtain analytical solutions. In this paper we 
also attempt to obtain the particular nonlinear solution 
for growing structure. Obviously, to proceed, additional 
assumptions are necessary. The problem becomes much 
simplified when wc adopt that the tangential velocity com- 
ponent is small relative to the radial one i.e. {ve/vrY <C 1. 
This results in a separation of eq. (14). The above assump- 
tion is consistent with the small value of the initial vortic- 
ity and the resulting equations reduce to the equation set 
(12-13) and (15), describing the rotationless flow of the 
magnetized fluid. 

To obtain a convenient form, we choose first the new 
time coordinate x = a~^/^ = t~^^^, instead of the cosmic 
time t. Our analysis becomes more clear when wc use the 
rescaled variables. Replacing the old functions 5, Wr, V' by 
the new variables X, Z2,Z4, where X = S, Z2 = arvr and 
Z4 = r^a^^ and defining the new independent variable 
^ = r/x, the continuity, the Euler and Poisson equations 
read 



4m 



1 



1 



6 



e+z2{o]x'{o = o, 



0, 



^2(0j4(0+^4(0 

2z4(0 + 54(0 



(16) 

(17) 
(18) 



where prime means the differentiation with respect to ^ 
and (3 = 4:TTG{poa^) (hereafter we take its value as /3 = |). 
This scaling is not only a mathematical manipulation. 
Introducing the ^ — variable reduces the differential com- 
plexity of the problem, allowing us to substitute the par- 
tial by the ordinary differential equations. In effect, the 
equations (16-18) govern the self-similar behaviour of the 
fluid motion, leading to the nonlinear 3'^'^-order equation 
for the rescaled velocity Z2{^) 



v(3) 



-15^1(0+4^^2(0 



tO 



z'm 

-24(0 
4(0 



le+z2{o 



ziiO = 0. 



(19) 



The detailed analysis of the above dynamical set 
(16-18) may be performed either through the qualitative 
methods of its asymptotic form or requires numerical in- 
tegration. 



Nevertheless, some general relation may be derived 
even before the solutions are obtained. The continuity 
equation (16) implies the generalized expression for den- 
sity contrast and velocity function Z2- 



ln(l + 



4(0 



e+z2{^) 



dO 



-dt. 



(20) 



This is a nonlinear counterpart of the known linear 
relation i.e. divv cx X (e.g. Willick & Strauss 1998), cou- 
pling p^fsv and thus B vs v. Lacking the analytical solu- 
tions we note that the asymptotic form for large velocities 
\z2\ ^ predicts the asymptotic property of the density 
contrast X + I \z2\~'^ ■ 

The linearization of the eqs. (16-18) and its result- 
ing solution for Z2 allows us to compare the time de- 
pendent part of Vr with the linear regime expressions 
v\^2 oc a^/^,o~^ (Peebles 1980). Finally, one has 



(4 



linx 



-e, (4'")2 



(21) 



and the implied linear solutions for density contrast are 

(X"")i = 9C^ (X"")2 = 4r^ (22) 

The analysis of the asymptotic behaviour of eqs. (16- 
18) at infinity (^ — > oo) becomes straightforward when 
we define the new variable C, ^ l/C- The subsequent 
linearization of the above equation set gives for ^ ^ cx) 
the following solutions for the velocity measure 



(4") 



-C" 



(4")2 = -C', (4")3 = const. 



and the density contrast 

(x"")i = (x"")2 = AC^. 



(23) 



(24) 



The above enables us to state that both velocity field 
and density contrast have solutions vanishing for ^ oo. 
Below we deal with the physical solutions i.e. regular at 
^ = and at ^ — > oo. 

4. Linear and nonlinear magnetic amplification 

The adiabatic matter compression determines the field 
amplification through the reached value of p. The amplifi- 
cation degree is controlled by the compression symmetry 
(see eq. (10)). For cylindrically collapsing structures one 
obtains the magnetic field amplification Wb = ^a^, thus 

taking into account eqs. (11) and (20) one has = 
const, while for example, in the case of spherical symme- 
try this quantity would be given by ^-^_|^^2/3 = const. On 
the other hand, the frozen field amplification expressed 
by the solution of the induction equation is defined by the 
velocity field given by eq. (19). To represent it directly by 
the time evolution of the velocity field we first present the 
velocity divergence in old variables (a;, r) 



V-V: 



1 d{rvr) 
r dr 



dz2 



dZ2 

dx 



(25) 
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Note that the 'self-similar' variable ^ is used here only 
for differentiation, since we hi 
exponent in eq. (11) becomes 

/ 



V-v , 1 
a 



2dz2 

X -——ax. 
ox 



£-2- The 
^ dx- 



(26) 



Replacing, in the above, the time variable by the scale fac- 
tor a and assuming the separability of the velocity function 
Vr = —V{r)v{a) (for collapse we have Vr{r,t) < ), one 
obtains the alternative form of eq. (26). 



dv{a) 
da 



+ v{a)a~^]da. (27) 



Thus the final expression for density and magnetic field 
compression due to nonlinear structure formation reads 



We 



X + 1 = exp 



V{r) 



I 



idv(a) . X _i 
a 2 ^ ' +v(a)a ^ 
oa 



da . 
(28) 



Combining the above equation, for the growing mode 
(f(a) oc a^/^) of the velocity field and eq. (21) yields the 
following amplification rate in the linear regime 



(29) 



The successive discussion of nonlinear behaviour will 
also deal with the growing mode. Knowledge of the ex- 
act nonlinear amplification Wb = exp[— / V-(-)d<] = 

J Ti 





'3V{r)a 




3 4 ■ 




■3 ; 


exp 


2r 


= exp 


-ra 


= exp 


2^ 



exp[ 



{2_|_22(j)d^] needs the solution of the very com- 
plex, quasilinear equation (19) for Z2- The zero value of 
the coefficient at the z^^ term allows us to determine 
the singularity points of this equation i.e. ^ = and 
+ Z2 =0. Since we postulate the regular velocity solu- 
tion at the structure centre we focus here on a later one, 
which appears when the ovcrdensity grows in the nonlin- 
ear regime. The singularity point position depends on the 
chosen boundary condition (a movable singularity). This 
is specific to the nonlinear differential equations (Bender, 
& Orszag 1978). Physically, the emergence of a singu- 
larity in the velocity and density field is an artefact of 
the collisionless fluid. This feature has also been found in 
A''-body simulations and the test-particle approach (e.g. 
Bertschingcr, 1985) to gravitational collapse, as an outer 
caustic which forms around the structure. (For coUisional 
matter, p ^ 0, it becomes a shock surface.) 

Since its analytical solution is difficult to obtain, be- 
low we present the results of numerical integration below 
the caustic surface. Independently we verify the asymp- 
totic behaviour of the solution at ^ ^ cxd. As expected 
from eq. (23), the velocity function -22 (^) and the density 
contrast X{^) approach zero at infinity. The principal fea- 
tures of the particular solution inside the structure and its 
asymptotic behaviour at infinity are illustrated in Fig. 1 
and Fig. 2 below. These figures show the density contrast 
and the velocity field as a function of ^ (or C) in the log-log 
scale. Both quantities are expressed respectively, in units 
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Fig. 1. The nonlinear excess of the density contrast X/X''" be- 
low the discontinuity (left panel), as a function of the ^-variable 
and (right panel) , above the discontinuity as a function of the 
C-variable. 
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Fig. 2. The nonlinear excess of the velocity measure ^2/^2° be- 
low the discontinuity (left panel), as a function of the ^-variable 
and (right panel) , above the discontinuity as a function of the 
C-variable. 



of (X"")i and (z^'")i or in units of (Xi'")2 and (z^'")2 at 
infinity. The presented profiles are only qualitative, in the 
sense that they are given by the numerical integration for 
a required set of the particular boundary conditions. The 
boundary conditions for numerical integration were taken 
from the linear solutions of eqs. (16-18). They clearly dis- 
play two distinctive regions corresponding to linear and 
nonlinear formation rates. Roughly, the linearity is repre- 
sented by the straight line (-^Un-)^ — 1 and the nonlinear 
region of the magnetic field and density contrast ampli- 
fication for ^ > 0.2 shows a strongly growing curve. It 
should be noted that the validity of the linear approach 
is slightly different in the case of density contrast and the 
velocity function^. 

The direct comparison of the linear and nonlinear 
amplification rate for a given structure radius shows 
that the nonlinear compression implies the magnetic field 
growth rate of several orders of magnitude. The predicted 
(cf. Fig. 1) amplification of the order of Wb ~ 10^ within 
the range (^1,^2) = (0.28, 0.48) corresponds to the evo- 
lutionary period given by |^ = (|y)^ ~ 6.5. Thus, the 
magnetic growth of the order of 2 x 10'^ may be attributed 
to the early nonlinear structure evolution between, for 
instance, redshift ~ 6 and redshift ^ 0. Certainly, the 
magnetic amplification becomes much greater as nonlin- 
ear evolution proceeds but then the validity of the applied 
assumptions becomes questionable. 



The density growth 



is faster than the velocity 
growth (jS^- It is clearly seen when the respective values 



(Xl.")i 



are compared i.e. 



(4'") 



-(0.2) = 1.01, ^ 



(0.2) = 1.02, 



(X' 



5^^(0.47) = 6.15 and 



■(0.47) = 2.18. 
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5. Conclusions 

The goal of this paper is to present the amplification rate 
of the magnetic field associated with the forming gravita- 
tional structure of cyhndrical symmetry. The widespread 
conviction that the large-scale structures are filled with 
microgauss cosmological magnetic fields motivates our in- 
terest in amplification processes during their evolution. 
On the other hand, we see a high degree of filamentarity 
in the galaxy redshift surveys (e.g. Sathyaprakash et al. 
1998). This demonstrates that we deal with magnetized, 
elongated structures of axial symmetry. 

To determine their magnetic structure growth, sev- 
eral simplifications are needed. We used here two cate- 
gories of simplifying assumptions: physical (i.e. p ~ and 
Fl ~ 0), constraining the results to the early nonlinear 
phase and geometrical ones — requiring the radial motions 
and thus the axial fields. On the basis of a formal solution 
of the induction equation we obtained the exact analyti- 
cal expression for linear field amplification. The relation- 
ship between the density contrast and the magnetic field 
strength is established through the velocity field diver- 
gence. However the major conclusion concerns the nonlin- 
ear phase. The magnetic field may be effectively enhanced 
there. The density contrast growth is stronger than the 
velocity, achieving the nonlinearity regime earlier. The ra- 
dial structure of the magnetic field and density contrast 
are identical, in general — nonhomogeneous. 

Contrary to this highly idealized model, in the realistic 
situation, the centrifugal forces will stop the collapse. This 
will however occur in the successive, virialization phase, 
when the matter will become collisional and then shocked. 
Therefore, a proper description requires more elaborate 
application of the fluid model. Within its current limita- 
tions the above applied symmetries seem to be less weighty 
than the physical assumptions. According to previous pa- 
pers (e.g. Siemieniec & Woszczyna 2004, Bruni et al. 2003) 
the more degenerate, pancake geometry leads to compara- 
ble amplification results. Introducing cylindrical symme- 
try enables us instead to depict the magnetic field profile 
inside the structure. The substantial enhancement of the 
matter density accreting onto collapsing structure indi- 
cates that significant magnetic fields may be produced in 
its outer region — the precursor of the future shock. 

Acknowledgements 

We thank the referee, Dr. G. Sigl for criticisms and 
useful comments. This work was supported in part by 
'Komitet Badah Naukowych' projects: 'Cosmic magnetic 
fields' (G. S-0) and 'Rigorous description of scalar and 
tensor perturbations in flat and open Priedmann uni- 
verses' (Z.A. G). 

References 

Barrow, J., Ferreira, P., & Silk, J., 1997, PRL, 78, 3610 

Bender, CM., & Orszag, S.A. 1978, Advanced Mathemat- 
ical Methods for Scientists and Engineers (McGraw-Hill 



Book Company, New York) 

Bertschinger, E., 1985, ApJS, 58, 39 

Bruni, M., Maartens, R., & Tsagas, C, 2003, MNRAS, 

338, 785 

Groth, E., & Peebles, P., 1975, A&A, 41, 143 

Kim, E., Olinto, A., & Rosncr, R., 1996, ApJ, 468, 28 

Lesch, H., & Chiba, M., 1995, A&A, 297, 305 

Lifshitz, E., 1946, J. Phys. (USSR), 10, 116 

Peebles, P., 1980, The Large-Scale Structure of the 

Universe (Princeton Univ. Press, Princeton) 

Peebles, P., 1993, Principles of Physical Cosmology 

(Princeton Univ. Press, Princeton) 

Sathyaprakash, B., Sahni., & Shandarin, S., 1998, ApJ, 
508, 551 

Siemieniec, G., & Woszczyna, A., 2004, A&A, 414, 1 
Wasserman, I., 1978, ApJ, 224, 337 
Willick, J., & Strauss, M., 1998, ApJ, 507, 64 
Zeldovich, Y., 1970, A&A, 5, 84 

Zeldovich, Y., Ruzmaikin, A., & Sokolov, D., 1980, 
Magnetic Fields in Astrophysics (McGraw-Hill, NY) 



